Diffractive arrays of gold nanoparticles near an interface: critical role of the substrate 
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The optical properties of periodic arrays of plasmonic nanoantennas are strongly affected by 
coherent multiple scattering in the plane of the array, which leads to sharp spectral resonances 
in both transmission and reflection when the wavelength is commensurate with the period. We 
demonstrate that the presence of a substrate (i.e., an asymmetric refractive-index environment) can 
inhibit long-range coupling between the particles and suppress lattice resonances, in agreement with 
recent experimental results. We find the substrate-to-superstrate index contrast and the distance 
between the array and the interface to be critical parameters determining the strength of diffractive 
coupling. Our rigorous electromagnetic simulations are well reproduced by a simple analytical 
model. These findings are important in the design of periodic structures and in the assessment of 
their optical resonances for potential use in sensing and other photonic technologies. 
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I. INTRODUCTION 



Scattering of light by periodic metallic structures has 
been well studied for over a century in the context of 
diffraction gratings P^ The subject has recently received 
renewed attention^ w ith the prediction^ and exper- 
imental observation^"^ f interesting optical phenom- 
ena that result from the interaction between the geo- 
metrical resonance associated with light diffraction and 
the excitation of localized surface-plasmon resonances 
in metallic nanoparticles, which play the role of plas- 
monic nanoantennas P 3 ^ In addition to the interesting 
physics revealed in such systems, a number of applica- 
tions hav e bee n prop osed, including nanoscale energy 
transport) 14 ! 15 ! sensing PSC3 an d modifying spontaneous 
emissionpO which rely on the improved quality factor re- 
sulting from the reduction in radiative damping of the 
array as compared to localized plasmons excited in iso- 
lated particles. Recent advances in the control of the 
angular emission from quantum dots are also based on 
diffractive coupling of antenna elements.^ 

Two-dimensional arrays of nanoantennas can be pro- 
duced with good fabrication contr ol by techniques such as 
electron-beam lithographyP"lA12Ql contact printing, 12 and 
colloidal chemistry!^ These samples are most commonly 
manufactured on a substrate of high refractive index 
compared to the upper medium (typically air, or water 
for biosensing applications). It has been suggested from 
experimental data^and theoretical modeling 22 that such 
an asymmetric configuration is incompatible with the ex- 
istence of delocalized surface modes, and therefore pre- 
vents the observation of efficient narrowing of the local- 
ized surface-plasmon lineshapes. This is in contrast to 
other studies of fluorescence in particle arrays P^ 

Here, we elucidate the conditions under which strong 
diffractive coupling may occur in the asymmetric config- 
uration. In particular, we examine the reflectivity of gold 



nanoantenna arrays in close proximity to a substrate. 
A sharp transition is observed when varying the parti- 
cle size above a certain threshold, leading to the sudden 
emergence of lattice resonances, and establishing a clear 
difference in the behavior of lithographically patterned 
nanoantenna arrays depending on the thickness of the 
metal layer. Likewise, we find the refractive-index con- 
trast and the array-substrate separation to be critical 
parameters. 



II. MODELING PARTICLE ARRAYS NEAR A 
SUBSTRATE 



The system under study is schematically depicted in 
Fig. [ij,. It consists of an infinite square array of gold 
spheres placed in water close to a glass substrate and illu- 
minated under normal incidence. For the particle dimen- 
sions under consideration, the optical response of each 
individual sphere shows a characteristic dipole plasmon 
resonance, as shown in Fig. lb for two different host ma- 
terials. The interaction between the particles in the array 
leads to additional, sharper features that correlate with 
the diffracted orders of the lattice, as shown for some 
specific geometrical parameters in Fig. lc. In particular, 
for a homogeneous environment, the <0,1> diffraction 
feature occurs when the wavelength in the medium is 
close to the period, X/n = a. This condition is obviously 
dependent on the refractive index n (Fig. lc, blue and 
green curves). Strikingly, the diffraction feature disap- 
pears when the array is placed in an asymmetric envi- 
ronment, close to an interface (Fig. lc, red curve). In- 
terestingly, the reflectance maximum associated with the 
dipolar plasmon shifts to the blue in this configuration as 
a result of interference between light scattered from the 
particles and the wave reflected at the planar interface 
(the specularly reflected beam at the planar interface in- 



terferes with the beam scattered from the spheres, which 
undergoes large phase shifts as the wavelength sweeps 
across the plasmon resonance, thus distorting the line- 
shape, unlike what happens in non-specular reflectance, 
for which such interference is not present). 

We set out to explain the origin of the suppression of 
diffraction close to the interface. The results of Fig. lc 
are obtained from rigorous multiple-scattering numeri- 
cal solution of Maxwell's equations,^ with the gold de- 
scribed by a tabulated measured dielectric functional and 
the rest of the materials having wavelength-independent 
refractive index. However, this method of solution does 
not provide insight into the origin of the effect, which we 
explore next using a simple analytical model that allows 
us to determine the most relevant parameters involved 
in the phenomenon of inhibition of collective resonances 
in the arrays. A variety of semi-analytical and numeri- 
cal techniques have been developed in this context G5H221 
For the geometrical parameters under consideration, it 
is reasonable to represent the particles as induced point 
dipoles, for which we obtain the polarizability a from the 
first Mie coefficient. - 3 

In the array, the dipole moment p-? at each particle site 
j satisfies the self-consistent coupled-dipoles equation 
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where E oxt is the external electric field and the sum gives 
the field induced by other particles j' at the lattice posi- 
tion Rj. Here, <?°(r) = (kf +VV) exp(ikir)/r is a Green 
tensor that yields the electric field produced by a dipole 
in the surrounding homogeneous medium of refractive- 
index rii an d light wave vector fci = kn\, where k is the 
free-space wave vector. 

When the wavelength in the medium 2-7T fk\ is com- 
mensurate with the period of the array a, scattering of 
the incident light by the array produces diffracted beams. 
Only some of these beams are propagating, but the re- 
maining evanescent waves play an important role. For 
wavelengths above the <0,1> threshold in the spectra 
of Fig. lc, only the specular beam is propagating, but 
four additional beams < ±1,±1 > become propagat- 
ing at shorter wavelengths. These beams are evanescent 
above the offset wavelength for diffraction, but in the 
presence of a substrate they may undergo multiple re- 
flections in the cavity formed between the array and the 
interface, thereby contributing to the reflectivity of the 
combined system. Thus our approach is to model the 
inhomogeneous environment using a Green tensor based 
upon the homogeneous environment and to add the effect 
of the substrate through (multiple) reflections from the 
substrate. 

The influence of the substrate on the beams is further 
illustrated in Fig. 2a, which presents dispersion diagrams 
in the form of reflectance as a function of the total and 
parallel light wave vectors, k\ and fen , respectively. The 
curves in Figure lc correspond to cuts of the color plots 



of Figure 2a along the vertical axis. For homogeneous 
environments (left and central plots), diffraction features 
are clearly visible near k x = lit ja. However, these fea- 
tures disappear in the presence of a substrate (right plot). 
Figure 2b shows that the <0,1> feature for normal in- 
cidence can receive contributions from diffracted beams 
with fen = 2n/a via momentum exchange with the lat- 
tice. Therefore, these diffraction beams have to be in- 
corporated into our model, and we show below that the 
substrate enhances their effect to the point of suppressing 
diffraction. 

It should be noted that diffraction originates in long- 
distance interaction between particles, which are polar- 
ized along the direction of the external field, x. Now, the 
dipolar field dies off faster than 1/r with the distance 
along x, and therefore, the < ±1,0 > beams are not ex- 
pected to play a significant role. Thus, in our minimal 
model we only include three beams with parallel wave 
vectors go = (0, 0) (specular beam) and g±i = ±(2-7r/a)y 
(diffracted beams). 

For an incident TE-polarized beam of parallel wave 
vector g m with unit electric field along x, the induced 
dipoles have the form p J = p m exp(ig m • Rj)x, which 
upon insertion into Eq. (IT]) yields 



1 



Pn 



1/a - G,, 



(2) 



where G m = J2j^o$° (Rj)exp(-ig m • R,-) is a lat- 
tice sum representing the collective interaction between 
dipoles (notice that the R = site is excluded from the 
sum).' 54 Both p m = p and G m = G are independent of m 
for the three g m beams under consideration. We show 



G = ya (Rj) 






(3) 



in the Appendix. 

The reflection coefficients of the array can be obtained 
upon examination of the scattered field. In particular, 
the x component reduces to 
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Now it is useful to expand the spherical waves of this 
expression in parallel wave vector space as 
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where k z = \/k\ — Q 2 is the normal wave vector compo- 
nent. Inserting Eq. (Is]) into Eq. Q and using the relation 
J2, exp(zQ • It,-) = (4n 2 /a 2 ) £j(Q _ g ), we find 
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FIG. 1: (a) Schematic representation of the system under study. An infinite periodic two-dimensional square array of gold 
nanoparticles is situated in a semi-infinite, homogeneous medium of refractive index n\ at a distance d above a substrate of 
refractive index n-z. The center-to-center separation between particles is a. Light is incident from the superstrate medium with 
a direction normal to the plane of the array and polarization along a lattice vector (x) . (b) Elastic scattering cross-section of 
an individual gold sphere immersed in water (n = 1.33) and in glass (n = 1.46), as obtained from Mie theory, (c) Same as 
(b) for a particle array of period a = 500 nm rather than an individual particle, also including the asymmetric configuration in 
which the spheres are in water with their surfaces located at a distance of 1 nm from a glass substrate. Vertical dashed lines 
indicate the position of the <0,1> and <1,1> diffraction conditions in the different media. All gold nanospheres are taken to 
have a radius of 35 nm. 



where g runs over reciprocal lattice vectors and k± = g± 
k z z. Noticing that g x = for the three beams included 
in our model, the scattered field can be approximated by 
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in which both k-t and k z depend on m through the lat- 
tice vector g m . The reflection coefficient of the array r a 
can be expressed as a 3 x 3 matrix with coefficients r"/ 
relating the incident beam f to reflected beams j. Using 
Eqs. Q and J6j, we find 



a 2 k z j 1/a — G' 



(7) 



where k Z j = \/k 2 — (27rm/a) 2 . 

The incident plane wave and the beams diffracted by 
the array are specularly reflected at the interface with 
the substrate, for which the coefficients of the reflection 
matrix r s are obtained from Fresnel's formula for TE 
polarization: 
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where k' z j — \Jk\ — (2irm/a) 2 . Multi-layered substrates 
can be straight-forwardly included in this analysis by in- 
troducing a suitable reflection coefficient. In particular, 



the case of substrates supporting surface modes may in- 
troduce new spectral features^; this is however beyond 
the scope of our present study. 

The composite array-substrate system forms an optical 
cavity, for which the total reflectivity must incorporate 
the effect of multiple internal reflections. A Fabry-Perot- 
type of analysis yields the combined reflection matrix 



r = r a -\ t a , 

1 j,aT> |.sp 1 



(8) 



where Pjf — Sjj> exp [ik Z jd] describes plane- wave prop- 
agation across the cavity and the transmission matrix 
of the array satisfies t%> — 5jj> + r^, . Finally, Eq. (TsT) 
gives an analytical expression for the specular reflectivity 
coefficient of the zero-order beam, 



'oo c 



,,2ikd 



[2e 2lk -- 



id a 



r 10 r 6l 



Ze r 01 r ll 



2 r 00 1J r 00 



e l r 00 r 00 



(9) 



This expression reveals an intricate dependence of the re- 
flectance of the composite system on the various physical 
parameters. The phase factors exp(2ik z jd) yield a pe- 
riodic modulation in regions far from intrinsic or lattice 
resonances (Fabry-Perot effect). The poles in Eq. (p)]) do 
not trivially expose the position of the modes supported 
by the structure because the numerator of Eq. Q may 
compensate for a possible divergence. Lattice resonances 



(a) 



water (n=1. 33) glass (n=1. 46) water/ glass 




k 2nla k 2n/a k 2n/o 



1 

10' 

io 2 
10 3 

<io- 4 



(b) 



;©v 



2n/a 

■HS* — 



-*£>' 



(c) 



-2n/a 



2n/a 



■+K 



6.10 



4.10 



0) 
U 

c 

-t-< 
U 



^2.10V,._ 




400 



500 600 700 

Wavelength (nm) 



800 



FIG. 2: (a) Dispersion diagrams showing the reflectance of 
square arrays of 35 nm-radius gold particles immersed in wa- 
ter (left), in glass (center), and in water at a surface-to- 
surface distance of 1 nm from a water-glass interface (right, 
see Fig. la). Here, fci is the wave vector in the host medium, 
feii is the parallel wave vector along a principal lattice direc- 
tion, and a = 500 nm is the lattice period. The incident po- 
larization is TE. (b) Scheme showing the first Brillouin zone 
(shaded area) and its vicinity. The dashed lines (light line and 
diffracted light line) correspond to the condition of grazing 
low-order diffracted beams. The circles signal the three beams 
contributing to the <0,1> diffraction spectral anomaly under 
normal incidence via (0, ±2tt /a) wave-vector transfers from 
the lattice (double-arrowed horizontal lines), (c) Calculated 
normal-incidence specular reflectance of the square array in 
the asymmetric environment considered in (a). Rigorous nu- 
merical solutions of Maxwell's equations including multipoles 
with orbital angular momentum up to I = 4 are compared 
with a simplified solution using only dipoles (Z = 1). These 
results, which require a large number of diffracted beams ~ 20 
to get convergence, are compared with our analytical model 
for dipoles [Eq. (pi] using only one beam (zero-order, red cir- 
cle in (b)) or three beams (all circles in (b)). 



may occur through two different sources. First, the coef- 
ficients r a contain a factor of the form (1/a— G) _1 , which 
is responsible for the lattice resonance in a self-standing 
array near the diffraction edge. Second, the reflectiv- 



ity coefficient r^ for the grazing diffractive orders has 
a factor k^ 1 = 1/ Jk\ — ^, which diverges right at the 

diffraction edge. It is precisely the interplay between the 
divergent terms in the numerator and the denominator 
of Eq. (pi) that is responsible for the cancellation of the 
diffractive coupling in an asymmetric configuration. 

This is clearly illustrated in Fig. 2c, in which we as- 
sess the accuracy of this model in describing the opti- 
cal properties of an array of small gold spheres near a 
substrate. The reflectance predicted by Eq. ^ (broken 
curves) is compared to a rigorous calculation based upon 
a multiple-scattering formalism described elsewhere^ 
(solid curves). The results from the analytical model 
are nearly identical to those of the full calculation near 
the <0,1> diffraction edge, and both predict a feature- 
less spectrum in that region. In contrast, a simplified 
version of the model accounting for only the specular 
beam [dotted curve; Eq. ffiih reduces to a scalar equation 
for the zero-order beam] shows a pronounced <0,1> fea- 
ture. The conclusion is clear: diffracted beams cancel the 
long-range coupling between particles in the presence of a 
substrate. Physically, this cancellation occurs because of 
the interference between the reflected grazing diffracted 
orders and the direct dipolar coupling in the superstrate 
medium. 

The accuracy of the three-beams model described by 
Eq. (p| is good near the <0,1> diffraction edge. How- 
ever, at shorter wavelengths, and in particular around the 
<1 , 1> diffraction edge, both the zero-order beam (dotted 
curve) and three-beams model (dashed curve) predict a 
non-existent spectral peak. This feature disappears when 
more beams are incorporated into the analytical model. 

Incidentally, we have also checked the validity of the 
dipolar approximation by including higher-order multi- 
poles in the rigorous calculation (up to I = 4). The re- 
sults are nearly indistinguishable regardless the number 
of multipoles present in the full calculation, as expected 
for spheres that are much smaller than the wavelength. 

In what follows, we investigate the influence of sev- 
eral physical parameters on the strength of the diffractive 
peak using the analytical model of Eq. (p| . 



III. RESULTS AND DISCUSSION 

In Fig. [3] we investigate the dependence of the re- 
flectance on the distance between the array and the sub- 
strate. The figure clearly shows the transition between 
the absence of diffraction when the particles are touching 
the interface and the sharp diffraction peak obtained at 
large distances, corresponding to the limit of a homoge- 
neous environment. As the array-substrate distance is 
reduced, this mode decreases in strength and begins to 
red shift for separations below 200 nm. Eventually, the 
mode disappears and a new, weaker peak is observed at 
the diffraction condition for the substrate medium. 

A similar transition occurs when the index of the sub- 
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FIG. 3: Calculated normal-incidence specular reflectance of a 
square array of 60-nm-radius gold spheres immersed in water 
and placed above a glass substrate (see Fig. la). The lattice 
period is 500 nm. The distance from the bottom of the spheres 
to the substrate is varied from 1 nm to 2 /im (see labels). 
Dashed curves: analytical model of Eq. Q. Solid curves: 
fully converged multiple-scattering numerical solution. 



strate is made to match that of the superstrate. In 
Fig. 4a we show reflectance spectra for a range of in- 
dex differences. As the index asymmetry is reduced, the 
mode strength increases dramatically after the asymme- 
try is below a threshold value. This is clearly observed 
in the evolution of the peak maximum with index con- 
trast (Fig. 4b, c). The threshold for diffraction inhibition 
seems to depend on both the period of the array and the 
size of the particles. In general, this threshold occurs at 
larger contrast when the particles are bigger or the period 
smaller, and the transition is smoother for smaller period. 
Actually, particles of larger size relative to the period 
deviate more from the ideal situation of small particles 
situated close to the interface, thus involving significant 
phase factors P [see Eq. ph] that make the diffraction fea- 
ture more robust against index contrast. A less intuitive 
result is that the array period does not substantially af- 
fect the maximum reflectance for a symmetric medium, 
although the fractional occupancy of the spheres is in- 
versely proportional to the period squared. This mode is 
associated with diffraction, and the reflectivity of the sys- 
tem is therefore strongly affected by coherent multiple- 
scattering, so that a simple geometric scaling rule is no 
longer applicable. 

Because spheres of increasing radius need to be placed 
further away from the substrate, the transition of Fig. [4] 
might perhaps be attributed to the effect discussed in 
Fig. [3] (the changing distance between the substrate 



and the sphere centers as the particle radius increases). 
To further discriminate between the two effects — the 
particle-centers separation from the substrate and the 
polarizability of the particles) — , we explore next arrays 
formed by elongated ellipsoids of increasing aspect ratio, 
showing that an increase in polarizability drives a tran- 
sition between suppression and emergence of diffraction 
features. Within the framework of the dipolar approxi- 
mation used in our analytical model, we adopt the po- 
larizability prescription developed by Kuwata-Gonokami 
et a/.p3 which provides a good approximation for sub- 
wavelength particles of moderate aspect ratio. Fig. [5ji 
presents results of this model for prolate gold ellipsoids 
immersed in water. The incident field is polarized along 
the long-axis of the particles. We have performed rig- 
orous T-matrix calculations 3 ^ to verify the accuracy of 
the model for the particles considered here. Elongat- 
ing the long axis of the particles enables us to tune the 
polarizability without altering the distance between the 
particles and the substrate. The scattering cross-section 
exhibits a localized surface plasmon resonance that is red- 
shifted and stronger as the aspect ratio increases. 

Using the polarizability of these particles in our analyt- 
ical model, we obtain the results shown in Fig. [5b for an 
array supported by a glass substrate. A transition very 
similar to that of Fig. [4^, is observed as the polarizability 
of the particles increases. The diffractive coupling in the 
superstrate medium is suppressed by the presence of the 
substrate for the 50 nm spheres, but as the long-axis of 
the particles is increased a new peak gains strength at 
wavelengths above the diffraction condition in the sub- 
strate. 

These results are qualitatively similar to those ob- 
tained for spheres of increasing radius. Fig. [6] represents 
the transmittance, reflectance, and absorbance of gold- 
sphere arrays immersed in water and supported on glass 
for various particle sizes and fixed pitch of the array. The 
vertical dashed line represents the onset of the <0 , 1> or- 
der diffraction. The arrays show an absorption and reflec- 
tion feature for a wavelength slightly to the right of the 
onset when the particle radius exceeds a value ~ 60 nm. 
This is the result of a lattice resonance involving the col- 
lective interaction of the spheres close to the condition for 
which the noted order of diffraction becomes grazing. 37 
The resonance is increasingly broadened and redshifted 
as the particle size increases. This effect is clearly visible 
when the radius is larger than ~ 80 nm. 



IV. CONCLUSION 

Our results clearly show that the presence of a sub- 
strate can reduce or even suppress diffraction in parti- 
cle arrays. This effect is important in experiments in- 
volving particles of small height relative to the period™ 
Diffraction can be recovered in samples with larger metal 
particles P21 These conclusions are important in the design 
of complex antennas involving interaction between metal 
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FIG. 4: (a) Calculated normal-incidence specular reflectance of a square array of 60nm-radius gold spheres and period equal 
to 500 nm immersed in a fluid and placed at a surface-to-surface distance of 1 nm from a glass substrate (see Fig. la). The 
refractive index of the fluid is varied in the range n\ = 1.2 — 1.5 around the glass index ri2 = 1.46. (b,c) Maximum of the 
reflectance spectra as a function of m for two values of the array period and three different sphere radii (see text insets). All 
spectra are calculated with the analytical model of Eq. q9p . 



parts at distances of several wavelengths on a substratePS 
They may also offer ways of performing sensing by de- 
tecting small variations of index of refraction in a fluid 
environment P^l 
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Appendix A: Lattice sum for normal incidence 
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We represent in Fig. [7] the calculated lattice sum G, 
as defined in Eq. (pf. This involves a poorly convergent 
sum that is accelerated by separating it in parts that are 
computed respectively in momentum and in real space, 
following the methods introduced by|39l The sum shows 
characteristic divergences at values of the wavelength for 
which a diffracted beam in the Rayleigh construction be- 
comes grazing, as explained elsewhere P^ 
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are indicated as vertical dashed lines. 
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FIG. 7: Real and imaginary parts of the lattice sum G as 
a function of wavelength A normalized to the period a for 
normal incidence on a square array. The wavelength is given 
relative to the medium in which the lattice is immersed. 



